Jamming as a critical phenomenon: a field theory of zero-temperature grain packings 
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A field theory of frictionless grain packings in two dimensions is shown to exhibit a zero- 
temperature critical point at a non-zero value of the packing fraction. The zero-temperature con- 
straint of force-balance plays a crucial role in determining the nature of the transition. Two order 
parameters, < z >, the deviation of the average number of contacts from the isostatic value and, 
<4>>, the average magnitude of the force per contact, characterize the transition from the jammed 
(high packing fraction) to the unjammed (low packing fraction state). The critical point has a 
mixed character with the order parameters showing a jump discontinuity but with fluctuations of 
the contact force diverging. At the critical point, the distribution of (f> shows the characteristic 
plateau observed in static granular piles. The theory makes falsifiable predictions about the spatial 
fluctuations of the contact forces. Implications for finite temperature dynamics and generalizations 
to frictional packings and higher dimensions are discussed. 

PACS numbers: 



CO 



I 

C 

O 

o 



m 

O 
in 
o 



-a 
c 

o 
o 



X 



Introduction In a remarkably diverse range of sys- 
tems, the transition from a flowing, liquid state to a 
jammed, solid state is heralded by a dramatic slowing 
down of relaxations^. Does an equilibrium critical point 
underlie this glassy dynamics ? The debate surrounding 
this question has been spurred by the absence of any ob- 
vious static signature accompanying the rapid increase 
of time scales Purely dynamical scenarios have been 
proposed^ to explain time-scale divergences with no ac- 
companying static divergences. For thermal systems, a 
different perspective has been offered within the frame- 
work of an avoided critical point 3] and a scaling the- 
ory based on the existence of a zero-temperature critical 
point In a more recent development, it has been sug- 
gested that the mechanism of jamming in both thermal 
and athermal systems is controlled by a zero-temperature 
critical point (J-point)||. 

Experiments on weakly sheared granular media indi- 
cate that at a critical packing fraction there is a transi- 
tion accompanied by slow dynamics, vanishing of mean 
stress, increasing stress fluctuations and a change in the 
distribution of contact forces 0,0- Simulations indicate 
a critical point occurring at zero temperature and a pack- 
ing fraction close to the random close packing value 
At this critical point the grain packing is isostatic, having 
reached the special coordination where all contact forces 
are completely determined by the packing geometry 
A theory based on this observation predicts a diverging 
length scale associated with the mechanical stability of 
the network 8J . In a different theoretical approach, an 
analogy has been drawn between the jamming transition 
and k-core percolation||. 

In the present work, it is shown that a statistical field 
theory of two-dimensional, zero-temperature, frictionless 
grain packings exhibits a critical point. This point sep- 
arates a disordered phase from an "ordered" one char- 
acterized by two order parameters: (i) the magnitude of 



the force per contact, <4>>, and (ii) <z>, the deviation 
of the contact number per grain from its isostatic value. 
At the critical point, the fluctuations around < 4>> di- 
verge but those around < z > go to zero. An analytic 
prediction for P(F), the distribution of contact forces, 
is in excellent agreement with experiments 0, and the 
theory makes falsifiable predictions regarding the spatial 
correlations of the forces. 

Statistical Ensemble In granular matter, which is 
athermal, a natural control parameter is the packing 
fraction[lJ. For short-range repulsive potentials, the cor- 
responding statistical ensemble is one with a fixed av- 
erage pressure [loj. The probability P[{i"j}] of a grain 
packing with the set of positions {r^} can be obtained 
by using the maximum entropy principle ^lj, which 
also forms the basis of the Edwards ensemble of gran- 
ular packings [l^. Maximizing the entropy, S[P] = 
~ S{ ri } -P[{ r i}] m -f[{ r i}] subject to the constraint of 
fixed average pressure leads to[l3j|: 



P[{r i }] = (l/Z)exp(ap({v i })) 



(1) 



where 2>({r;}) = (l/^)X)y r ij W^a is the pressure of 
the configuration {r^}( U is the interaction potential). 
The "canonical" partition function, Z(a) — ^2' r -P[{ r i}} 
is the generating function of all statistical averages. The 
Lagrange multiplier a plays the role of inverse tempera- 
ture: a = —dS/d<p> and the prime on the summation 
restricts it to grain configurations satisfying the equa- 
tions of mechanical equilibrium for frictionless packings 
in d dimensions: 



dN eqs : ^ F t , 







<z> N/2 eqs : F tj = /(r„) 



(2) 



(3) 



Here < z > is the average number of contacts per grain, 
Fij is the magnitude of the contact force between grains 
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the contact forces through: Fj 



hj/ where j' and 



FIG. 1: Illustration of height field. The dashed lines 
are contact forces and the heights associated with the 
loops are represented by the arrows encircling a loop 



i and j, and /(r.jj) is the function specifying the inter- 
grain force law. At the isostatic point, <z>— z iso = 2d, 
the number of equations in Eq. 0is exactly equal to the 



number of unknown forces, 



and, therefore, 



the forces are uniquely determined by these equations. 
The positions {r{\ can then be obtained by inverting Eq. 
13 For <z>> Zi SO , Eqs. [3 and are coupled and the 
coupling can be parametrized by e which encodes how 
sensitive the contact forces are to changes in the grain 
positions [l5|: 
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(4) 



For hard spheres, e = 0, the two equations are decoupled 
and the only packings for which jF^-j's and {ry}'s can be 
determined are the isostatic ones|l4|. For e << 1, small 
variations in grain positions can lead to force changes 
comparable to the average force and Eqs. [21 and |3| are 
weakly coupled. An effective way of calculating the par- 
tition function, in this limit, is to sum over all the so- 
lutions, {Fij}, to Eq. [3 for a given set {r,} and impose 
a Gaussian constraint; e _( - e / 2 ^^ Pi ~ Pi - > . Here pi is the 
pressure calculated from {ri} and p { ( is calculated from 

{Fay- 



Pi 

Z(a) 



Pi = 
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e («i: l p.) e ~^/2i: l (p.-p°) 2 _ ( 5 ) 



The sum over Fij is constrained to only those forces 
which satisfy force-balance. These constraints can be in- 
corporated through the introduction of a set of auxiliary 
variables, the loop forces [lfij|. Loop forces or 'heights', 
{hi}j0l are vectors associated with the the voids en- 
closed by grains. The height vectors are related to the 



are the voids bracketing the contact ij (cf Fig. 0. 
Since the Fy-'s around a grain sum to zero, the map- 
ping of forces to heights is one to one, up to an arbi- 
trary choice of a single height. For frictionless grains, 
the requirement of Fy being parallel to ry leads to re- 
strictions on the heights. In developing the field theory, 
the heights are coarse grained over a mesoscopic region 
with a length scale much larger than the grain radius 
but smaller than a typical length scale over which the 
height fields var y Ha to define a continuous height field. 
It can be shown[17) that the constraint on the coarse- 
grained field resulting from the frictionless property is 
that the height field be divergence free and, in 2d, this 
implies that the height field can be written in terms of 
a scalar potential[16| ip- h x — d y ip ; h y = —d x tp. A 
similar approach led to the definition of a potential for 
frictional grain-packings [lq . The stress tensor, coarse 
grained over a mesoscopic region of area A, around the 
point r: o"°(r} = (l/A)J2 jcA J2k r jkFjk, is expressed in 
terms of ip as|l6| 



-d x d y ip 



-d x d y ip 



and the coarse grained pressure is, p (r) = Tra — 
V 2- i/>(r). The partition function can, therefore, be written 
in terms of the unconstrained field ip 

The coarse-grained field ip carries with it a weight 0,[ip] 
which counts the number of microscopic {F^} configura- 
tions giving rise to the same ip field. Arguments similar 
to the ones employed in height-maps of loop models[l8|. 
lead to Q[ip] - e ~(i/2) j d d r(v 2 i,)\ The basic reason _ 

ing is that configurations which are completely flat, i.e. 
V 2 i/> = 0, have the largest number of loops (dashed in 
Fig. ^) along which the contact forces can be shuffled 
without violating the force-balance constraint and, there- 
fore, they are favored entropically. 

At the isostatic point, z = Zi SO , p = V 2 ?/'- For packings 
close to the isostatic one, an expansion in z — Zi SO leads 
to 



p(r) = V 2 ijj + (z-z, 



SO A 

Az 



V 2 V> + (z — z iso ) 



v 2 v> v 2 v> 



(6) 



The second set of equations follow by noting that the 
extra pressure, ^ \z Uo , due to the introduction of an ad- 
ditional contact at the isostatic point is the pressure at 
this packing divided by the number of contacts. Just 
as the mapping from Fij — > ip, the mapping from ri to 
the coarse-grained pressure, p(r), involves the calcula- 
tion of a weight u>[p\. It can be argued that u>[p] ~ 
exp(— K\(S7p)\ 2 ) by picturing the packing of N grains 
in a volume V. Entropically, it is favorable to have the 
grains occupy as much of the volume as possible, lead- 
ing to small spatial variations of the packing density and 



3 



p(r). Large pressure gradients occur in the entropically 
unfavorable packings where there are large variations in 
the packing density. 

With these weights, the partition function is: 



<z> 



Z 



c/2f d d r( P (r)-p°(r))f 
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3 0q4^qi- 



(qi +Q3) ' (Q2 +Q4)] 



4 qi^q2 ( /^q3 ( /'q4 c 'qi+q2+q3+q4J (7) 

Here, cf> = V 2 ip represents the magnitude of the force per 
contact (c/Eq. EJ, and the field z has been redefined to 
z — Zi S o which is restricted to the set of integers. The pa- 
rametes a, e and K have been scaled to absorb resulting 
factors of Zi SO , 

Critical Point The Hamiltonian H reflects the com- 
petition between favoring = and the "field" a 
favoring non zero </>. It provides a model for studying 
the response functions of grain packings with small but 
finite e. We assume that the integer restriction on z can 
be ignored ^ as long as e ^ 

In investigating whether or not there is a finite-a phase 
transition involving the vanishing of one or more order 
parameters, the fields in the Hamiltonian in Eq. [7] are 
expanded around their average values; q =<</> q > +Cq, 
Zq =< Zq > +?7 q and the hamiltonian written as H = 
H (< cj) q >, < Zg >) + Hi(< <j> q >, < z q >;( q> r)q). The 
order parameters < 4> q > and < z q > are obtained by 
minimizing the effective potential, T(< 4> q >,< z q >) = 
#o(< g >, < z q >) - ln(jrn q dC q n q d7? q e" ffl ). The sim- 
plest, non-trivial approximation, is obtained by calcu- 
lating the fluctuations, < |C q | 2 > and < |^ q | 2 >, at 
the loop level, replacing all 4-point averages by 2-point 
aver ages |20j. and assuming spatially uniform order pa- 
rameters: < 4> >=< 4> q =o >, < z >=< z q= Q >. To leading 
order in e: 

< ICql 2 > _1 = l + e<z> 2 -1/ <(f)> 2 +K<z> 2 q 2 

(8) 



< Kl" > _i - £<cf» 2 +K <(f» 2 q 2 

Incorporating these results in the effective potential and 
minimizing with respect to the order parameters leads 
to: 



<(j)> -a(z iso + <z>) + e <z> 2 <4» + 



1 



<0> 







e(<0> 2 + 



1 - 1/ <cj)> 



r ) <z> -a <4>>= (9) 



Solving these equations near a c = 2/z isoi the order pa- 
rameters behave as: 

«f» = (a/a c )(l + (l-a c /a) 1/2 ) 



(a/e)(l-a c /a) 1 / 2 



(10) 



For a > a c , there is an ordered phase characterized by 
two order parameters. For a < a c , T(<0>,<z>) ceases 
to have any local minima or maxima (cf Fig. [3Jl and 
<4>> jumps discontinuosly to 0: the physical limit of its 
allowed values. From Eqs. ISland HUI asa^a c , <fi—^l and 
the q — force fluctuations diverge: < cf) 2 > — < <j>> 2 ^ 
(1 — a c / a)^ 1 / 2 . This type of transition is indicative of 
the end of a line of metastable equilibrium similar to 
spinodal critical points. Unlike spinodals, however, the 
transition is accompanied by the disappearance of any 
local minimum; a phenomenon observed in models with 
rigid constraints such as certain dimer models 21]. The 
fluctuations of z vanish as < z 2 > — < z > 2 ~ (1— a c /a) 1 / 2 . 

The structure factor of the magnitude of the contact 
forces (the <j> field), 5(q|) =< |C q | 2 >, provides detailed, 
measurable information about the spatial fluctuations of 
the contact force. Eq. |H1 predicts 



S(q) 



1 



t/a{e,K) + q 2 t 2 



where f ~ (1 - 1/ < <t» 2 ) ~ (1 - a 2 /a 2 ) 1 ' 2 and a(e, K) 
is a characteristic length scale of the model. The length 
scale £ vanishes at the critical point and S(q — 0) diverges 
as l/£. There is another length scale related to the width 
of the structure factor peak: Z -1 = go where go = 1/Vo?- 
The width diverges as (1 — a 2 /a 2 )^ 1 ^. The appearance 
of these two length scales, and the specific form of S(q), 
are predictions of the theory that should be testable in ex- 
periments and simulations. The length scales, I and £, are 
associated with correlations of contact forces and their 
vanishing implies that the length scale over which the 
grain packing behaves as an elastic solid is going to zero 
or, equivalently, the length scale over which the packing 
is floppy is diverging. The two exponents, 1/2 and 1/4, 
associated with the divergence of these complementary 
length scales are reminiscent of the exponents discussed 
in theories and simulations of J-point jjj, 0, • 

Distribution of contact forces Experiments 0, Q and 
simulations E| have shown that changes in the dis- 
tribution of P(F) are associated with transitions involv- 
ing the vanishing of stress. In the field theory, < <j> > 
corresponds to F and the distribution, P(< <j) >) ~ 
e -n«f»,<z(<4»)). The critical regime has < z(< <f> > 
)) >~ {cf Eq. U implying T(< (f> >) ~ (1/2) < > 2 
— az iso < cj)> +ln(< <f) >)■ The ln(< <f> >) term results 
from integrating out the z field and embodies the phys- 
ical effect of large contact number fluctuations for small 
contact forces. This term leads to a power-law decay 
P(< (f>>) ~ 1/ < <f) > near </> = and provides a qual- 
itative explanation of the experimental dataQ. Fig. [21 
clearly demonstrates that the approach to a c is accom- 
panied by changes in the low-force regime of P(F) with a 
peak giving away to a plateau as the force fluctuations di- 
verge. Since diverging force fluctuations imply vanishing 
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FIG. 2: The distribution P(<cf>>) of contact forces (from 
bottom to top) a — a c — 0.7, 0.5 and 0.025. The inset 
shows the <z>= cut of the potential, T(<cf)>,<z>) 
at different a with a — ■+ a c from bottom to top. 

of the shear modulus ~ l/(<0 2 > — <<j)> 2 ) , the shape 
change is connected to the disappearance of yield stress. 
The theory does not yield the exponential tail of the force 
distribution which is not unexpected in a coarse-grained 
model since the large force behavior results from a stress 
redistribution at the particle level 0, ■ 

Conclusions A field theory, founded on general prop- 
erties of T = packings of frictionless grains in 
two-dimensions, exhibits a critical point separating a 
jammed, "ordered" phase with finite yield stress at large 
packing fractions (a > a c ) from an unjammed phase 
where the order parameters vanish. A hallmark of the 
transition is a change in the distribution of contact forces 
reflecting diverging force fluctuations. Two length scales, 
associated with the spatial fluctuations of the contact 
forces, are predicted to go to zero at the critical point 
and should be detectable in experimental measurements 
of the structure factor of contact forces. 

If the zero-temperature critical point controls the low 
temperature and weakly-driven behavior, then glassy dy- 
namics follows from general scaling arguments \A l23| . At 
finite temperatures, the force-balance constraint is vio- 
lated and the height-field has defects 01- The interac- 
tions between these defects determine the stability of the 
zero-temperature critical point and, therefore, it will be 
crucial to understand these interactions. The extension 
of the theory to frictional packings is within reach since 
the loop- force formalism exists[l^|. Although the simi- 
larities between two and three dimensions, observed in 
simulations, suggest that an extension to higher dimen- 
sions is possible, this remains an open question. 
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